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With the commissioning of new powerful facilities with high production rates the experi- 
mental investigation of eta decays will become a precision tool for testing low-energy QCD 
in the SU (3) flavour sector. The anticipated numbers of 10 s - 10 9 observed etas per year 
at DA$NE, ITEP, CELSIUS and GRAAL will allow for a test of chiral perturbation the- 
ory (ChPT, see Weinberg (1979) and Gasser and Leutwyler (1984, 1985)) predictions up 
to next-to-leading order or even beyond. Experiments will include high-accuracy measure- 
ments of the main decay modes r\ — > it-kit and rj — > 77, as well as the detection of rare 
decay modes predicted by ChPT. The mode 77 — > 77 involves the anomalous Wess-Zumino- 
Witten (WZW, Wess and Zumino (1971), Witten (1983)) action and has been thoroughly 
analysed in combination with the analogous tt° and 77' decays in theory as well as in exper- 
iment. However, the WZW action also generates more complicated interactions, which are 
of particular interest for the following reason: at 0(p A ) the WZW lagrangian completely 
covers the odd intrinsic parity sector of ChPT and is the only piece in the 0{p 4 ) ChPT 
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lagrangian which has actually predictive power by itself, because its low-energy coefficient 
is completely determined by the structure of the WZW functional and, in contrast to the 
Gasser-Leutwyler lagrangian for the even intrinsic parity sector, need not be determined 
by means of phenomenological input. Hence, the validity of the predictions derived from 
the WZW action is of fundamental interest for the validity of ChPT and QCD itself. The 
process 77 — > 7r + 7r~77 has the remarkable feature that it involves three different types of 
WZW interactions in a consistent 0(p i ) treatment according to Weinberg's power counting 
(Weinberg (1979)). Such a calculation is a tree- level calculation, because an odd number of 
pseudoscalar mesons participates in the process and, thus, one vertex in the Feynman dia- 
grams contributing to the process must have odd intrinsic parity. The invariant amplitude 
for the processes 77 — > 7r + 7r _ 77 and 77 — > 7r°7r°77 is constructed from the standard and 
the anomalous C$zw lagrangian, 
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where only those terms of C w \ w have been retained which give a contribution to the two pro- 
cesses under investigation. The covariant derivative is defined as D^U = d^U + ieA fl [Q, U], 
where Q = diag (2, — 1, — 1) /3 is the diagonal quark charge matrix. The current quark 
masses are contained in M = diag(m u , m d , m s ), and the constant B is related to the quark 
condensate. Since we want to include 77-77' mixing at a phenomenological level, we assume 
nonet symmetry. Hence, the matrix U contains the pseudoscalar octet fields as well as the 
eta singlet field. We use the exponential representation of the pseudoscalar meson fields, 
U = exp(i$/F$), where F$ is the pseudoscalar meson decay constant and the field matrix 
is defined as 
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The three types of Feynman diagrams for the charged decay mode are a diagram with a four- 
meson interaction from & 2 \ a tt° or i] propagator and a one-meson-two-photon vertex from 
Cwzw (cl ass 1 diagram), a diagram with a three-meson-one-photon interaction from C^) zw , 
where another photon is emitted from a charged pion line (class 2 diagram) and a contact 
diagram with a three-meson-two-photon interaction from C^zw The class 1 diagrams are 
gauge invariant by themselves, whereas for class 2 and 3 only the sum of the diagrams 
is gauge invariant. It is interesting to note that the matrix element for r\ — > 7r + 7r~77 is 
closely related to that of 77 — > 7r + 7r~7r° (see Bos (1994) for a treatment in ChPT). Thus, 
a measurement of the decay could allow for a consistency check of the process with two 
photons in the initial state. The invariant amplitude for the neutral decay mode receives its 
only contributions from class 1 diagrams. 

Calculating the photon spin sum for the invariant matrix element squared and integrat- 
ing over four of the five independent variables of the four-particle phase space results in the 
diphoton spectrum dT/dz (z = s 7 /m^), where s 7 is the invariant cm. energy squared of 
the diphoton system. More details concerning the calculation can be found in Bijnens et 
al. (1992) and Knochlein et al. (1995). The diphoton spectrum of the charged decay mode 
(Fig.l (left)) clearly reflects a bremsstrahlung singularity for s 7 — > 0. Another singularity of 
the invariant amplitude due to the diagram with the propagating 7r° influences the spectrum 
in a very narrow energy region around s 7 = m^ and could be removed by introducing the 
final width of the 7r° in the propagator. In the rest of the spectrum the influence of the dia- 
grams with the propagating neutrals is almost invisible. The partial width T(i] — > 7r + 7r~77) 
is calculated as a function of an energy cut 5rrib rems applied around s 7 = 0. Our result for 
Snibrems = 195 MeV is smaller than the experimental upper limits (Price et al. (1967)) by 
about three orders of magnitude. However, for Snibrems ~ 30 MeV our calculation yields 



a branching fraction of approximately 1CT 4 which should be within the reach of the new 
facilities. The diphoton spectrum of the neutral decay mode (Fig. 1 (right)) is clearly 
dominated by the 7r° pole diagram which is present due to isospin and G parity violation 
at the four-meson vertex. The eta pole diagram, which is the only diagram present in 
the limit of conserved isospin, plays a minor role in the dynamics of the decay. The partial 
width for the neutral mode, which is obtained by integration of the diphoton spectrum and 
application of a reasonable energy cut around ^J~s^ = m^o, is by more than two orders of 
magnitude smaller than for the charged mode. Hence, it will probably be extremely difficult 
to detect this mode. 
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FIG. 1. Left figure: Diphoton energy spectrum dT/dz (z = s 1 /mi) for the decay 
7] — ► 7r + 7r~77. The dotted line is the full calculation, the dashed line is calculated without class 1 
diagrams and the solid line is proportional to the phase space integral. 

Right figure: Diphoton spectrum for the decay r] — > 7r°7r°77. The dash-dotted line is the prediction 
of chiral perturbation theory for m u = 5MeV and = 9MeV. The dotted line is the ChPT 
prediction for m u = rrid = 7MeV. The dashed line is a calculation, where we determined the 
strength of the ?77r 7r 7r interaction from the decay rj — > 7r°7r°7r and of the ijij'^tt interaction 
from the decay rj — ► ry7r°7r°, including also an additional r/g^s 71 " 071 " interaction. The solid line is 
proportional to the phase space integral. 



